

JU 

- /. , ,» - A'VS 


AIAA 2000-0757 

Wake Vortex Transport and Decay in Ground 
Effect: Vortex Linking with the Ground 

Fred H. Proctor and David W. Hamilton 
NASA Langley Research Center 
Hampton, VA 

Jongil Han 

North Carolina State University 
Raleigh, NC 


Z *-!# 1 • 

2D-I’ ^ 

ID- 


fmz 

. «■ , 


r * • ' A. #» 

¥ . . wF « * W 


38th Aerospace Sciences 
Meeting & Exhibit 

January 10-13, 2000/ Reno, NV 


For permission to copy or republish, contact the American Institute of Aeronautics and Astronautics 
1801 Alexander Bell Drive, Suite 500, Reston, Virginia 20191-4344 



AIAA-2000-0757 


Wake Vortex Transport and Decay in Ground Effect: Vortex Linking 

with the Ground 


Fred H. Proctor* and David W. Hamilton* 
NASA Langley Research Center 
Hampton , Virginia 


Jongil Han* 

North Carolina State University 
Raleigh, NC 


Abstract 

Numerical simulations are carried out with a three- 
dimensional Large-Eddy Simulation (LES) model to 
explore the sensitivity of vortex decay and transport in 
ground effect (IGE). The vortex decay rates are found to 
be strongly enhanced following maximum descent into 
ground effect. The nondimensional decay rate is found to 
be insensitive to the initial values of circulation, height, 
and vortex separation. The information gained from 
these simulations is used to construct a simple decay 
relationship. This relationship compares well with 
obsen'ed data from an IGE case study. Similarly, a 
relationship for lateral drift due to ground effect is 
constructed from the LES data. In the second part of this 
paper, vortex linking with the ground is investigated. Our 
numerical simulations of wake vortices for IGE show that 
a vortex may link with its image beneath the ground, if the 
intensity of the ambient turbulence is moderate to high. 
This linking with the ground (which is observed in real 
cases) gives the appearance of a vortex tube that bends to 
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become vertically oriented and which terminates at the 
ground. From the simulations conducted, the linking time 
for IGE appears to be similar to the linking time for 
vortices in the free atmosphere; i.e., a function of ambient 
turbulence intensity. 

Nomenclature 
B aircraft wing span 

b u initial vortex separation - k BI 4 

g acceleration due to gravity 

L Xr Ly domain length in x and y directions, respectively 
p, domain depth 

M mass of generating aircraft 

r radius from vortex center 

r r radius of peak tangential velocity 

t time coordinate 

T nondimensional time - t V 0 / b a 

T at Z = Zmirt 

y a airspeed of generating aircraft 

V G initial vortex descent velocity - TJ (2 n b 0 ) 
x horizontal coordinate along flight path 

y horizontal coordinate lateral to fight path 
Y nondimensional lateral coordinate - y/b n 

z vertical coordinate 

minimum altitude achieved during vortex descent 
Z, initial height above ground of wake vortex 
F vortex circulation 

r o initial vortex circulation - M g!(b 0 p V A ) 

p air density 

€ ambient turbulence (eddy) dissipation rate 

^-component of vorticity 
v kinematic viscosity 

rj nondimensional eddy dissipation - 
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I. Introduction 

Aircraft wake vortices are formed when vorticity 
generated by lift rolls up into two primary trailing 
vortices. The two vortices have a characteristic initial 
separation, b„, which is proportional to the wingspan of 
the generating aircraft. The wake vortices are considered 
to be “in ground effect” (IGE) when located within one 
initial separation of the ground (i.c., z < b 0 ).‘ The 
prediction for the location and intensity of these vortices 
is important for systems that are to manage safe aircraft 
spacings. 

The complex influence of the ground on wake vortex 
behavior has been the focus of numerous investigations, 
including the observations by Harvey and Perry, 2 * 
Halloed Kopp, 4 Hallock and Burnham, 5 Burhnam and 
Hallock, 6 and Rudis et al.J as well as the numerical 
studies by Schiling, 8 Robins and Delisi, 9 Zheng and Ash, 10 
Corjon and Poinsot, 11 Corjon and Stoessel, 12 and Proctor 
and Han. 13 Qualitative understanding of basic vortex 
transport mechanisms for IGE is known, but quantitative 
relationships are needed. Wake vortex decay is known to 
be enhanced by proximity to the ground, but much is to be 
learned about this behavior 6 

Three-dimensional instability of wake vortices near 
the ground is largely unexplored. Lateral linking of the 
vortex pair associated with Crow instability 14 should be 
less likely due to the increased separation forced by the 
ground’s presence. Ground linking -- the linking of a 
vortex with its image vortex beneath the ground plane -- 
has been observed, 15 but thought to be rare. 16 Linking 
instabilities can increase the uncertainty of a vortex’s 
position and may affect the vortex decay rate as well. 

Semi-empirical vortex prediction algorithms have 
been developed by Robins et al. 1 that include the ground’s 
influence on lateral separation and vortex rebound. 
Similar algorithms 17 have been incorporated within 
NASA’s Aircraft VOrtex Spacing System (A VOSS). 18 19,20 
Improved accuracy of these prediction algorithms depend 
on the better understanding of vortex transport and decay 
near the ground. Numerical experiments with a Large 
Eddy Simulation (LES) model are being conducted in 
order to provide guidance for the enhancement of these 
prediction algorithms. 21 

This paper is a continuation of the three-dimensional 
LES study described in Proctor and Han, 11 with further 
focus on vortex decay within IGE. In Proctor and Han, 
we investigated the sensitivity over a wide range of 
ambient turbulence levels of a wake vortex for a landing 
L-I011 that was observed on 26 September 1997 at 
Dallas-Fort Worth (DFW) airport 22 The simulations 


showed only a weak sensitivity to the level of ambient 
turbulence, and were in very-good agreement with Lidar- 
derived observations. The results indicated that enhanced 
decay from ground effect begins a few seconds after the 
vortices descend to their minimum height, with vortex 
decay prior to these times being primarily influenced by 
stratification and ambient turbulence. The study also 
found that vertical oscillations associated with vortex 
rebound were less prominent for higher levels of ambient 
turbulence, since the secondary vortices were more likely 
to be diminished by the stronger turbulence. Also 
addressed was the influence of ground stress on wake 
vortex transport. The study determined that viscous 
stress reduced the velocity near the ground and acted to 
decouple the primary vortices from their sub-ground 
images. If ground stress was ignored, unrealistic 
divergence of the vortex pairs ensued. The study, did not 
address Crow instability, or similar linking instabilities. 

Table 1. Salient Characteristics of TASS 


Primitive equation / non-Boussinesq equation set Time- 
dependent, nonhydrostatic, compressible. 

Meteorological framework with option for either three- 
dimensional or two-dimensional simulations. 

Large Eddy Simulation model with 1 st-order subgrid 
scale turbulence closure - Grid-scale turbulence 
explicitly computed, while effects of subgrid-scale 
turbulence modeled by Smagorinsky model with 
modifications for stratification and flow rotation. 

Ground stress based on Monin-Obukhov Similarity 
theory. 

Optional boundary conditions - open conditions utilizes 
mass-conservative, nonreflective radiation boundary 
scheme. 

Explicit numerical schemes, quadratic conservative, 
time-split compressible- accurate as well as highly 
efficient, and essentially free of numerical diffusion. 
Space derivatives computed on Arakawa C-grid 
staggered mesh with 4 th -order accuracy for convective 
terms. 

Prognostic equations for vapor and atmospheric water 
substance (e.g. cloud droplets, rain, snow, hail, ice 
crystals). Large set of microphysical-parameterization 
models. 

Model applicable to meso-yand microscale atmospheric 
phenomenon. Initialization modules for simulation of 
convective storms, microbursts, atmospheric boundary 
layers, turbulence, and aircraft wake vortices. 
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since the assumed domain had a limited axial depth. In 
the present paper, we address the sensitivity of vortex 
decay for IGE, by conducting LES of wake vortices with 
different initial heights, separations and circulations. The 
results are normalized and a simple analytical 
representation for IGE decay is given. Similarly, a simple 
model based on the LES results is proposed for lateral 
drift due to ground effect. In the second part of this paper, 
we examine the occurrence of ground linking and its 
sensitivity to ambient turbulence. 

II. The Model and Initial Conditions 


Thermodynamic Equation (Potential Temperature ): 
dO _ I d0p o uj ^ Q dp ()U j 

& Po Po d *j 


I 


Po*Xj 


3 r v ^0 , 

lPo K ti 3 — 1 


} Xj 


with the Potential Temperature defined as: 

p Kj 

e=r(-^)<; 


The numerical model used in this study is a three- 
dimensional LES model called the Terminal Area 
Simulation System 23 (TASS), which has been adapted for 
simulation of wake vortex interaction with the 
atmosphere. 24 The numerical model used in this study 
differs from that in many previous investigations, in that: 
1 ) time -dependent computations are carried out in three- 
dimensional space, 2) the formulation is essentially free of 
numerical diffusion, 25 3) the computations are LES and at 
high Reynolds number, 4) the experiments are initialized 
with realistic turbulence fields, and 5) realistic surface- 
stress and subgrid turbulence-closure formulations are 
assumed. Salient characteristics of TASS are listed in 
table 1. 

Model Equations 

The TASS model contains a prognostic equation set 
for momentum, temperature and pressure, and employs a 
compressible time-split formulation. Omitting moisture and 
coriolis terms (which are not used in these simulations), the 
TASS equation set in standard tensor notation is as follows: 


In the above equations, u, is the tensor component of 
velocity, t is time, p is deviation from atmospheric pressure 
P , r is atmospheric temperature, p is the air density, C p and 
C v are the specific heats of air at constant pressure and 
volume, g is the earth’s gravitational acceleration, R (i is the 
gas constant for dry air, P oa is a constant equivalent to 1 000 
millibars (10 5 pascals) of pressure. Environmental state 
variables, e.g., p„ , P„ and 0 o , are defined from the initial 
input sounding and are functions of height only. 


A modified Smagorinsky first-order closure is used for 
the subgrid eddy viscosity as: 


11 


du, / du, | °Uj 2 () in .2 
d X) d.xj <>x, '3 dxk 


Km = l s 

x sjl - a, Ri s - u 2 Ri r 


The subgrid eddy viscosity for momentum, K M is modified 
by the Richardson numbers, for stratification, Ri s , and for 
flow rotation, Ri r , with oq =3 and 0b= 1 .5. 


Momentum: 

d Ul | H dp 
dt Po d Xi 

1 d 

+ — 3 — P 0 Km ( 

Po d *j 


d -^ U -+m^ iL + g(H- 1)8,3 

OXi V.I 


Xj 

dm 


KKj 


Ulj 


d xj dx { 


2 

3 dx k 


Sul 


Buoyancy Term: 

0() PoCp 


Pressure Deviation: 

dp C„P du j 


dt Cv d 


Xj 


= p 0 g itjSji 


The subgrid length scale, l s , is determined from the 
grid volume and is matched to the appropriate length scale 
close to the ground where the flow is under-resolved. That 
is: 


/,= 


aA 

aAJJHaAJk zTjj 
l+( aA/k z C 

kz 


z>aA/k 
aA/k> z> Az/2 
z<Az/2 


where k is von Karman’s constant, and where m and a are 
invariant constants with values defined as m = 3 and a = 
0.15. The filter width is based on the minimal resolvable 
scale: 

A=[ 2A.v 2A_v 2Az ] in 
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where Ax, Ay, and Az are the numerical grid sizes in the 
respective x, y, z direction. 

The ground boundary is impermeable with nonslip 
velocity specifications. The surface stress due to the 
ground is determined locally from the wind speed, surface 
roughness, and the local thermal stratification. Details of 
the surface formulation are in the appendix of reference 

[13]. 

Turbulence Initialization 

Prior to vortex initialization, an initial field of 
resolved-scale turbulence is allowed to develop under an 
artificial external forcing at low wavenumbers. 26 The 
approach is similar to that in recent studies w-ith TASS, 
where wake vortex decay and the development of Crow- 
instability is examined within a Kolmogorov 27 spectrum 
of homogeneous turbulence (Han et al. 2i ' 29 ). The method, 
however, is slightly different for the present study, due to 
the inclusion of the ground. Since periodic boundary 
conditions are assumed only at the horizontal boundaries, 
the turbulence forcing is applied only to horizontal 
velocity over each horizontal plane. Nevertheless, the 
influence of the horizontal two-dimensional forcing 
spreads quickly to the vertical direction as well as to the 
vertical velocity through the mass continuity. 



Figure 1 . Turbulence energy spectrum at z= 14 m before 
vortex injection. Here, subscripts 2 and 3 denote 
crossflow and vertical directions, respectively. 


Because the TASS code uses a finite difference 
numerical scheme, the forcing is achieved by performing, 
first, a two-dimensional fast Fourier transform (FFT) at 
every large time step, then adding a constant amplitude to 
all the modes w ith integer wavenumbers w r hose magnitude 


is less than 3.0, and finally, performing an inverse FFT 
back to the physical space. At the same time, the 
horizontal domain average temperature and velocity fields 
are forced to maintain their initial vertical profiles by 
subtracting the difference every time step. Due to subgrid 
dissipation, the simulation can reach a statistically steady 
state, in the sense that the mean turbulence kinetic energy 
oscillates in time around a constant value. 

Figure l shows a one-dimensional energy spectrum 
with a -5/3 slope of Kolmogorov’s 27 spectrum when the 
turbulent flow Field has achieved a statistically steady 
state. Figure 1 indicates that our approach can produce a 
well-developed turbulent flow field that possesses 
Kolmogorov’s inertial subrange. The eddy dissipation rate 
E is estimated from the well-known technique of fitting 
Kolmogorov's theoretical spectrum in the inertial 
subrange to the simulated spectra. 

Vortex Initialization 

The initial wake vortex field is specified with a 
simple vortex system that is representative of the post roll- 
up, wake-vortex velocity field. The vortex system is 
initialized with the superposition of two counter-rotating 
vortices that have no initial variation in the axial direction. 
The tangential velocity, V, associated with each vortex, is 
determined from: 

/ s0.75 

l-Expt-10\^\ I 

w here, r is the radius from the center of the vortex, D* 
the circulation at r >> r r , and B is the span of the 
generating aircraft. The above formula is based on Lidar 
observations of wake vortices measured early in their 
evolution. The model is applied only at r > r c and is 
matched with a Lamb model profile for r < r r . The 
values assumed for initial vortex separation and 
circulation are derived from an aircraft’s weight, span, 
and airspeed according to conventional formula based on 
elliptical loading; i.e., b„ = nB/4, and 
Appropriate vortex image conditions are applied to the 
initial wake field to ensure consistency and mass 
continuity at the TASS model boundaries. 

Except when otherwise noted, the initial vortex 
parameters (table 2) are taken from the observed aircraft 
parameters for the L-1011, as was used in Proctor and 
Han. 13 

All simulations are conducted with dimensional 
variables and assume turbulent flow with a rotational 
Reynolds number (Tyv) of - 10 7 . 
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Table 2. Initial vortex 

parameters (baseline). 

Parameter 

Value 


Vortex spacing (b 0 ) 

37 m 


Generating height (Zj) 

16 m 

(0.432 b a ) 

Vortex circulation (T 0 ) 

390 m 2 s 

1 

Vortex core radius 

3 m 

(bj 12.3) 


Model Domain Parameters 

Two domain sizes are used for the experiments in 
this paper. A short domain, which has been truncated in 
the A-direction (L x , L y , L z = 81 m, 370 m, 81 m), and a 
long domain (L Xt Ly , L z - 451 m, 337 m, 81 m). For an 
assumed initial vortex separation of b„=37 m, this would 
translate into: (2.2 x 10 x 2.2) b 0 for the short domain and 
(12.2 x 9.1 x 2.2) b„ for the long domain. Both long and 
short domains are resolved by the grid sizes listed in table 
3. 

Table 3. Grid Resolution. 

Parameter __ Value 

Vertical resolution (Ac) 1.5 m 

Lateral resolution (Ay) 1.5 m 

Axial resolution (Ax) 2.0 m 


The advantage of the short domain is that it allows 
economical calculations of vortex behavior and simplifies 
the analysis of results. The longer domain, however, is 
less likely to hinder three-dimensional linking instabilities. 
A comparison of results from the two domain sizes is 
shown in the appendix. 

Determination of Vortex Position and Circulation 

Based on the vortex position in each crossflow plain, 
circulation is computed for each crossflow plane according 
to: 

r r <*)= l! r Cx dyfo 


The circulation is not computed for any plane where the 
vorticity vector is beyond 30" of the x-axis. 


A 5-/5 meter averaged-circulation, is computed 
according to: 


r u> = 


SlTrdr 

1 > 


where a =5 m , and b = 15 m. The 5-15 m average 
circulation is chosen to characterize the intensity of the 
vortex, since it quantifies the hazard faced by an 
encountering aircraft 30,3 1 


A mean average circulation is reported for each time 
interval by computing a mean of the circulations from each 
crossflow plane: 


' N; ' 


a,b 


Similarly, a mean vortex position is reported for each 
time interval by averaging the positions from each 
crossflow plane. 


HI. IGE Sensitivity 

Three sets of experiments are conducted in order to 
determine the sensitivity of vortex decay w ithin IGE. The 
short domain is assumed which will inhibit the occurrence 
of linking instabilities. All of the simulations assume 
neutral stratification, no mean ambient winds, and an 
ambient turbulence dissipation rate of e =9.6 x 10 * m 2 
s' . Other initial conditions for the simulations are listed 
in tables 2 and 3. 


Sensitivity to Initialization Height 

Five experiments are conducted for a vortex 
initialization height (Z,) of: 0.32 b n , 0.5 b 0 , 0.65 b,„ 0.84 
b,„ and 1.0 b„. Comparisons of the experiments are shown 
in Fig. 2 as function of nondimensional time. Note that 
the time coordinate for each experiment is offset by the 
time of maximum descent into ground effect; i.e., by 
T ( ,~T(zmin )■ The curves for normalized 5-15 m average 
circulation, as well as the normalized lateral positions tend 
to collapse, once the offset time T (; is accounted for. As 
also shown in our previous study, 13 vortex decay is 
significantly enhanced following maximum penetration 
into IGE. In addition, these results show the decay rate to 
be more or less independent of the vortex initial height, 
with the enhanced level of decay beginning roughly 0.25 
nondimensional time units following 7>, 


A formula for the decay of average circulation can 
be obtained from these results as 

r -2(T -T ) 2/J 

L =Exp! i oo j 1 ( 1 ) 

1 oo 5 


where T o0 = T Cl + 0.25 (i.e., 0.25 nondimensional time 
units following the time of maximum penetration into 
ground effect), and L (to is the circulation at 
nondimensional time T oa . The above formula is 
independent of ambient turbulence and can only be 
applied for T > T (; +0.25. 

Differentiating the above gives the rate of decay as: 
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4 


L i L 

Too dT 


WT~T 00 ) 


1/ 


-Exp{ 


-2(T~T 00 ) 


2/3 


A plot of the curve from Eq. 1 is shown also in Fig. 
2c and agrees very well with the LES results. 

An empirical relationship for lateral spread within 
IGE can be determined from the LES results as well. The 
vortices, which start at a lateral distance from the flight 
path of y = ± ftbo , diverge lateral with time due to the 
influence of the ground. Figure 2b indicates (in the 
absence of crosswind) that the vortices asymptote toward 
y = 2b„\ approaching a maximum separation of about 4b n . 
Based on these results, an empirical relationship for 
vortex drift due to ground effect is: 

Y= 1.385 (T-Trf 227 for T > To + 0.25 (2) 

where Y is the normalized lateral position from the flight 
path. The vortex separation, b = 2Yb ot also can be 
predicted for IGE with Eq. 2. A plot of the curve from 
Eq. 2 is compared with LES data in Fig. 2b. 

The lateral drift rate can be determined by 
differentiating (2) giving: 

V G = — = 0.3173(T-T C f om for T > T r , + 0.25 
dT 

where V (i is the nondimensional velocity for lateral drift 
due to ground effect (this velocity can be made 
dimensional by multiplying by the initial vortex sink rate, 

V 0 ). 

Sensitivity to Initial Vortex Spacing 

In order to further evaluate the sensitivity of our 
results and test the validity of Eq. 1, three additional 
experiments are conducted, which assume different initial 
vortex spacings ( b „ ) as listed in table 4. An initial 
circulation of 400 m 2 s l and an initial height of Zj= 0.432 
b n is assumed for each experiment. Values for Z, and 
nondimensional turbulence dissipation are slightly 
different in these experiments due to the variation in b iv 


Table 4. Initial vortex separation ( bj and corresponding 
initial altitude (ZJ and dimensionless turbulence intensity 
frj) for sensitivity ’ to initial vortex spacing cases. 


b c (m) 

Zi (m) 

T1 

30 

13 

0.0671 

37 

16 

0.0897 

49 

21.2 

0.1291 


IGE Sensitivity Test for Varying Vortex Heights ( r}=0.091 ) 
Vortex Vertical Position History (Port) 

(A) 



IGE Sensitivity Test for Varying Vortex Heights ( n=0 091 ) 
Vortex Lateral Position History (Port) 



IGE Sensitivity Test for Varying Vortex Heights (n=0.091 ) 
Vortex 5-15 m Averaged Circulation History (Port) 



Figure 2. Sensitivity to vortex initiation height. 
Nondimensional height (a), nondimensional lateral 
position (b), and nondimensional 5-/5 m average 
circulation (c) v.v nondimensional time. “ Decay model ” in 
(C)from Eq.I and “ Curi e fit ” in (B) from Eq.2. 
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IGE Sensitivity Test for Varying Initial Separation (b 0 ) 
Vortex Vertical Position History ( Zi=0.43b 0 ) 

1.4 

L (A) 

1.2 

- 

"-a 1 

* 08 

|o.6 

< 


0.4 

0.2 

b g ■ 30 in, r] ■ 0.0671 

— b 0 - S7 m, r, • 0.0887 

"l ijiI. i ■ ii) • * i i i i i iiiiiii 

°C 

1 2 3 4 5 


Nondime nslonal Time (t b^VJ 

IGE Sensitivity Test for Varying Initial Separation (b 0 ) 
Vortex Lateral Position History { Zi=0.43b 0 ) 

2.5 

r (B) 

"“a 2 

f 


|15 

M 


* « 

Jr 

2 1 

a 

3 

0.5 

b„ * 30 m, n a 0.0671 

L b„ * 37 m. n - 0.08*7 

Modal FK 

' . . . » l . . . . 1 . . • i 1 i ■ ■ i » i * < » 1 

°c 

1 2 3 4 5 


Nondime nstonal Time (t b/VJ 

IGE Sensitivity Test for Varying Initial Separation <b 0 ) 
Vortex 5-15 m Average Circulation History ( Zi=0.43b 0 ) 

1.2 

r (°) 

1 

0.8 

••••*•••••• b„ m 30 m, n ■ 0.0S71 

V b„ » 37 m, n * 0.0887 

1 ; b„ * 49 m, r| a 0.1 291 

rn, Decay Modal 

^0.6 

u 


0.4 


0.2 

■ . . . . i . . . . i . . . . i . . - * i . * » » i 

°C 

> 1 2 3 4 5 


Nondime nstonal Time (t b^VJ 


Figure 3. Same as Fig. 2, but for sensitivity to initial 
vortex spacing. 


IGE Sensitivity Test for Varying Initial Circulation 



Figure 4. Same as Fig. 2, but for sensitivity to initial 
vortex circulation. 
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The results presented in Fig. 3 show that the curves 
again collapse to a similar value when normalized, and 
that the proposed formulas for IGE decay and lateral drift 
are in excellent agreement with the LES results. 
Enhanced levels of decay due to IGE begin at 7=0.5 
(which is Ta+0.25) for all three cases. 

Sensitivity to Initial Circulation 

In the third set of experiments. Five simulations with 
different initial circulations are conducted (table 5). All 
else is assumed equal, with r\ being different due to its 
dependency on T 0 . 

Again, the curves collapse with a tight spread (Fig. 
4), and the proposed IGE decay and drift formulas are in 
excellent agreement with LES results. 


Table 5. Initial circulation (Tj and corresponding 
dimensionless turbulence intensity for T 0 sensitivity cases. 

r o (mV') 

n 

250 

0.1420 

300 

0.1 183 

350 

0.1014 

400 

0.0887 

450 

0.0789 


IV. Comparison with Observed Case. 

In this section, the new formula for IGE decay 
(Eq. I ) is compared with the LES simulation for a landing 
L-I01J , which was observed at 20:09 UTC on 26 
September 1997 at Dallas-Fort Worth (DFW) airport. 
The input sounding for temperature and crosswind is 
shown in Fig. 5. The simulation is conducted with the 
short domain and with the parameters listed in tables 2 
and 3. The simulation includes an ambient turbulence 
field with, e =1.654 x 10 m 2 .v" ? (q =0.2349), which is 
very close to the ambient value observed at z~40 m. The 
average circulation from both the port and starboard 
vortices is show in Fig. 6. The slightly faster decay for 
the starboard vortex may be due to the opposite sign 
vorticity of the ambient crosswind. Also included in Fig. 
6 for comparison, is measured Lidar data for the port 
vortex. The LES results and Lidar data show reasonable 
agreement with the circulation decay given by Eq. 1 . 

V. Ground Linking Experiments 

In the free atmosphere, the linking process begins as 
the counter-rotating circulations of a vortex pair connect, 
producing crude vortex rings. The time at which this 
linking takes place is a known function of the turbulence 


intensity, or specifically, the eddy dissipation rate 28,32,33 
The goal of the experiments in this section is to examine 
ground linking and any sensitivity that it may have to the 
level of ambient turbulence. All of the experiments 
assume the tong domain so as to permit linking 
instabilities. The experiments also assume the parameters 
listed in tables 2 and 3, as well as the initial sounding in 
Fig. 5. The nondimensional turbulence strengths, q, used 
in these experiments, range from 0.233 to 0.75. These 
values represent a typical range between moderate to very 
turbulent atmospheric boundary layers. 



temperature and b) crosswind for DFW , 26 September 
1997, 20:09 UTC . 


DALLAS 1997 IN GROUND EFFECT CASE 
97/09/26 20:09:24 GMT 

Vortex 5-15 m Averaged Circulation Comparison 



Figure 6. Average circulation vs nondimensional time 
for obsen'ed DFW case. Comparison of LES, Lidar, and 
IGE decay model ( Eq. 1 ). 
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Results from these simulations showed obvious 
ground linking to occur when ri exceeded 0.3 (e g. Fig 
7). For lower values of nondimensional turbulence 
intensity, pronounced vertical oscillations developed 
without obvious linking between the wake vortex and 
its ground image (e.g. Fig. 8). In none of the cases 
(which are all initiated well within 1GE) did lateral 
linking occur. Influence of the very weak crosswind 
shear is unknown (see Fig. 5b), but ground linking only 
occurred with the port vortex. 

The visualization of the simulated wake vortex as 
it links with the ground is shown in Fig. 7. The port 
vortex becomes nearly vertical as it links with the 
ground just after T =2. For comparison, a photograph of 
an actual ground link is shown in Fig. 9. 



Figure 9. Observed case of ground linking (from 
NASA Langley photo archives, photograph L 90- 
02903 ) 


Both port and starboard vortices in Fig. 7 rapidly 
dissipate after linking, becoming undistinguishable 
from background turbulence at T=4. Note that in the 
vicinity of where the ground connection is to occur, the 
starboard vortex is transported upward and the port 
vortex downward. 

Figure 8 depicts the wake vortex system for a 
more moderate level of turbulence. Although an 
obvious link with the ground did not occur, a major 
portion of the port vortex remained parallel and very 
close to the ground. Again, note that the starboard 
vortex is at a relatively high altitude where the port 
vortex is closest to the ground. 


A ground linking factor can be defined similarly to 
Crow's linking parameter 14 as: 


I ) — £ min 

*-nuvc t’tnin 

where z mLX and z mm are respectively, the maximum and 
minimum altitude of one of the vortices. Either the port 
or starboard vortex may be considered linked with its 
ground image when the linking factor exceeds 0.85. 
Table 6 shows nondimensional times for p to exceed 
both 0.75 and 0.85 as a function of nondimensional 
turbulence intensity. Also included is the linking time 
as predicted by Sarpkaya’s 17 34 3 * re-derivation of Crow 
and Bate’s 32 theoretical formula for vortex lifespan in 
the free atmosphere. In our experiments, ground 
Unking did not occur for r| < 0.388 , although the 
vertical linking parameter did grow with time. 
Additional experiments were carried out for larger 
domain sizes assuming r\ -0.233, but little sensitivity 
was noted. Therefore, the linking time (or lack of 
linking) was not believed to be affected by our domain 
size. The influence of the crosswind shear on the 
suppression of linking for the lower turbulence levels is 
yet to be examined. 


Table 6. Sensitivity of Ground Linking to Ambient 
Turbulence (port vortex). 


T1 

Tb=0.75 

Tb=0.85 

1 Sarp 

0.233 

4.3 


2.4 

0.30 

2.90 


2.0 

0.388 

1.90 

2.10 

1.71 

0.50 

1.47 

1.63 

1.43 

0.75 

0.98 

1.14 

1.11 


Figure 10 shows a comparison of the time to 
ground link with experimental data 33 and theory 32 34 for 
vortex linking in the free atmosphere. The three 
simulations that had ground linking are in close 
agreement with Sarpkaya’s formula, 34 33 even though 
the theory w'as developed for lateral vortex linking 
away from any ground influence. 
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Vortex Lifespan (time to link or burst) 
vs Turbulence Intensity 



Figure 10. Linking time vs n and imetisional turbulence 
intensity. Time of ground linking from LES given by 
triangles, all other data for lateral linking (or vortex 
bursting ) out of ground effect (OGE). 

VI. Summary and Conclusions 

A study of wake vortices initialized in-ground 
effect has been performed using a validated three- 
dimensional LES model. The simulations, which 
included environmental turbulence, show strong decay 
of wake vortices following maximum penetration into 
1GE. Wake vortex linking with its ground image was 
simulated for moderate to strong values of atmospheric 
turbulence. 

The primary conclusions of this study are: 

1 . Wake vortex transport and decay for IGE can 
be nondimensionalized with conventional 
parameters. 

2. Appropriate normalization of parameters 
reduces wake vortex sensitivity to initial 
circulation, initial separation, and generating 
height. 

3. Vortex decay for IGE has minor sensitivity to 
the ambient turbulence level. 

4. Based on the LES results, a formula for IGE 
decay is proposed. This formula, with few 
dependent parameters, should be easy to 
incorporate in wake vortex prediction models. 

5. Also based on the LES results, a model for 
lateral drift rate due to ground effect is 
proposed. 

6. Both the decay and drift models show very 
good agreement with the LES sensitivity 
results. 


7. Ground linking is influenced by the level of 
ambient turbulence. Similar to lateral linking 
in the free atmosphere, the time to ground link 
is a function of the nondimensional eddy 
dissipation. However, ground linking did not 
occur for moderate to weak levels of 
turbulence. 

8. The increased lateral separation of vortices 
from ground effect can suppress lateral vortex 
linking. 

Future work needs to address the sensitivity of 
crosswind shear to the onset vortex linking, both in and 
out of ground effect. 
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Appendix - Sensitivity of Domain Size 

A number of experiments assume a truncated 
domain in the x-direction, so as to reduce computation 
time, inhibit crow instability and simplify analysis of 
the experiments. The consequence of shortening the 
domain is evaluated here. Two simulations are carried 
out with everything identical except for a short domain: 
(L Xt L y , L z = 2.2 b n , 10 b„, 2.2 b„f and a long domain 
(L x , Ly, L z ~ 1 2.2 b„, 9.1 b 0 , 2.2 b a ). The ambient 
turbulence level is assumed as r\=0.2349 and each are 
resolved by the grid sizes listed in Table 3. The 
ambient temperature and crosswind are from profiles 
observed at 20:09 UTC on 26 September 1997 at DFW 
(Fig. 5). 

Comparison of the average circulations from the 
two simulations with different domain sizes is shown in 
Fig. A-l. Lidar observations from the L-101 1 port 
vortex are also included. A comparison of lateral 
positions is shown in Fig. A-2. For this case, the 
figures indicate little sensitivity to domains size. 
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Figure A-l. Time evolution of the 5-15 m averaged 
circulation for port vortex from short domain 
( L x ~2.2b (> ) and long domain ( L x -12.2bJ simulations. 
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